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ABSTRACT 

We calculate the dissipation efficiency of relativistic reconfinement shocks. Building 
on previous work (Nalewajko & Sikora 2009), we consider different distributions of the 
external pressure. The average dissipation efficiency Cdiss is a function of the product 
of two parameters - the jet Lorentz factor Fj and the reconfinement angle 0r, which 
is related to the opening angle 0j and the external pressure index 77. The spatial dis- 
tribution of the dissipation rate strongly depends on 77. We discuss the significance 
of these results for the properties of relativistic jets in gamma-ray bursts and active 
galactic nuclei and propose that reconfinement shocks may explain a very high dis- 
sipation efficiency of the former and a moderate dissipation efficiency of the latter. 
Finally, we estimate the dissipation efficiency of the reconfinement shock associated 
with the quasi-stationary knot HST-1 in the jet of radio galaxy M87 and show that it 
is roughly consistent with the observational constraints. 
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1 INTRODUCTION 

Relativistic jets stand behind the brightest cosmic phenom- 
ena: gamma-ray bursts (GRBs) and blazars, a subclass of 
active galactic nuclei (AGNs). Their extreme isotropic lumi- 
nosities, up to ~ 10^^ erg s"^ for the former (GRB 080319B; 
Racusin et al. 2008) and up to ~ 10^° erg s~^ for the latter 
(3C 454.3; Abdo et al. 2011), cannot be plausibly explained 
without the relativistic Doppler effect. However, even taking 
this into account, in order for this radiation to be produced 
in the co-moving reference frame, a substantial fraction of 
the jet mechanical power, including the particle rest energy 
flux, needs to be dissipated, then transferred into a popula- 
tion of ultra-relativistic particles in a non-thermal accelera- 
tion process and finally radiated away through non-thermal 
radiative mechanisms. The total efficiency of these processes 
can be estimated observationally if the total jet power is 
known. In the case of GRBs, it can be well constrained by 
energetics of the afterglow phase and the radiative efficiency 
of the prompt phase has been claimed reach values up to 
~ 90% (Zhang et al. 2007). In the case of blazars, these esti- 
mates are less certain, but typical values for their luminous 
class of Flat Spectrum Radio Quasars (FSRQs) are ~ 10% 
(Celotti & Ghisellini 2008). Regardless of efficiencies of the 
particle acceleration and radiative processes, efficiency of the 
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energy dissipation in relativistic jets must be at least com- 
parable to these observational constraints. 

The most widely discussed means of energy dissipation 
in relativistic jets are shock waves, magnetic reconnection 
and instabilities. Shocks can arise within a jet when two 
regions propagating with substantially different bulk veloc- 
ities collide with each other. Such internal shocks provided 
the basic framework for theoretical models of blazars {e.g., 
Blandford & Konigl 1979) and GRBs {e.g., Rees & Meszaros 
1994). However, these models have been questioned on the 
grounds that they cannot account for required dissipation 
efficiencies. This is especially clear in the case of GRBs, 
for which several alternative models have been recently pro- 
posed, based on magnetic reconnection {e.g., McKinney &: 
Uzdensky 2010; Zhang & Yan 2011) or relativistic turbu- 
lence [e.g., Narayan & Kumar 2009). In the case of blazars, 
detailed calculations showed that a substantial contrast of 
initial Lorentz factors must be assumed (Spada et al. 2001). 
However, the occurrence of such a velocity contrast cannot 
be verified with current models of jet formation and accel- 
eration. Moreover, if the velocity modulations are related to 
processes at the black hole horizon scale, the internal shocks 
model predicts a particular length scale, a fraction of a par- 
sec, over which such shocks develop. There are now several 
arguments for the bulk emission of luminous blazars being 
produced at much larger distances from the central black 
hole {e.g., Sikora et al. 2008; Agudo et al. 2011). 
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The other possibility for the shocks is that they result 
from the interaction between the jet and its environment. In 
the case of GRBs, the jet is a transient phenomenon and has 
to plough through its host star and the interstellar medium, 
forming an external shock that dominates during the af- 
terglow phase {e.g., Meszaros & Rees 1997). In the case 
of blazars, the jet is relatively persistent and propagates 
roughly along a tunnel drilled over time, so perpendicular 
external shocks are not usually considered. However, the ex- 
ternal medium can exert a substantial pressure on the jet 
boundary, forcing it to recoUimate and triggering a recon- 
finement shock. 

Reconfinement shocks were first discussed by Sanders 
(1983) in the context of the kpc-scale jet of the radio galaxy 
NGC 315. The first analytical models were introduced by 
Canto et al. (1989) in the non-relativistic regime applicable 
to the jets of young stellar objects (YSOs) and by Komis- 
sarov & Falle (1997) in the relativistic regime. A unique 
signature of reconfinement shocks is that, unless the jet or 
external medium parameters vary significantly on the dy- 
namical time scale, they would be observed as stationary 
patterns. Daly & Marscher (1988) interpreted a stationary 
knot in the pc-scale jet of radio quasar 4C 39.25 as a nozzle of 
the reconfinement shock. More recently, a stationary^ knot 
HST-l has been discovered in the jet of radio galaxy M87 
at the 100 pc scale (Biretta et al. 1999) and subsequently 
it underwent a spectacular multiwavelength outburst {e.g. 
Harris et al. 2006). Stawarz et al. (2006) showed that the 
association of this feature with a reconfinement shock is con- 
sistent with both the properties of the host galaxy and the 
estimated jet power. However, short variability time scales 
required a very compact emitting region. Bromberg & Levin- 
son (2009) showed that efficient focusing of the shocked jet 
flow is possible, but requires substantial cooling of the post- 
shock plasma. Reconflnement shocks were also studied in 
the context of GRBs (Bromberg & Levinson 2007). 

The problem of dissipation efficiency of relativistic re- 
confinement shocks was first studied in Nalewajko & Sikora 
(2009), hereafter Paper I. It was found that the dissipation 
efficiency e^iss depends strongly on the product of the jet 
Lorentz factor Fj and the opening angle Gj. Here, we gener- 
alise this result, taking into account different distributions 
of the external pressure. We also show how this result can 
be applied to both GRBs and the jets of active galactic nu- 
clei. Because GRB jets are characterised by wide opening 
angles, reconfinement shocks provide a natural explanation 
of their high radiative efficiency in the prompt phase. In 
AGN jets, the efficiency of reconfinement shocks is much 
lower, because coUimation by a continuous medium limits 
the opening angle. We also estimate the efficiency of the re- 
confinement shock associated with the HST-l knot in the 
jet of M87 and show that it is roughly consistent with the 
observed luminosity of this radio galaxy. 

In Section 2, we present our simple model of the struc- 
ture of relativistic reconfinement shocks. The dependence of 
the dissipation efficiency on model parameters is discussed 
in Section 3. In Section 4, we discuss the applications of 




Figure 1. Geometric parameters of our reconfinement shock 
model. The jet is symmetric around the z axis and propagates 
from its origin at z = towards the reconfinement point at z = z^. 
The opening angle ©j , the reconfinement angle 0r and the max- 
imum jet width rm arc indicated. 



these results to astrophysical relativistic jets. Conclusions 
are given in Section 5. 



2 THE MODEL 

Reconfinement shocks result from the interaction between 
a jet and its surrounding medium. The simplest model of 
such a problem involves a cold, unmagnetized, spherically 
symmetric jet of Lorentz factor Fj, opening angle Qj and 
total power Lj ; and a static medium of pressure distribution 
given by Pc{z) oc , where < 2 is a constant and z is the 
coordinate measured along the jet axis. Figure 1 shows the 
geometric parameters of the reconfinement shock front of ra- 
dius rs{z), including reconfinement length Zr, reconfinement 
angle Or and maximum jet width Vm', as well as the contact 
discontinuity of radius rc{z). Parameters measured immedi- 
ately upstream and downstream of the reconfinement shock 
are denoted with subscripts 'j' and 's', respectively. 
The system of shock jump equations is 



j3s COs(6's — Qs) 

Usps sin(6ls - Qs) 



u^Ws sin^(6ls — as) + Ps 



FsUsWs sin( 



/3j COs{ej - Qe) , (1) 

UiPi sin(6lj - as) , (2) 

pjC^ sin {Oj - as) , (3) 

TjMjPjC^ sin(ej - Qs) , (4) 



where (3 = v/c is the dimensionless velocity, u = F/3 is the 
dimensionless four- velocity, w = pc^ +p + eis the relativistic 
enthalpy in the comoving frame, p is the mass density, e is 
the thermal energy density, as is the inclination of the shock 
front with respect to the jet axis and 0j,s are the inclinations 
of the velocity vectors. It is assumed that pj = 0. Given all 
the parameters of the upstream plasma, this system can be 
solved when the post-shock pressure ps is given. In Paper I, 
we noted that the structure of the shocked zone, the region 
between the shock front and the contact discontinuity, can 
be quite complex^. In principle, Ps{z) < Pc{z), so that the 
transverse pressure gradient can focus the post-shock flow. 
However, since the results on the dissipation efficiency pre- 
sented in Paper I are not very sensitive to the treatment of 
the shocked zone, we use the simple 'Model 1' from Paper 
I and assume that Ps{z) = Pc{z). The main improvement is 
that we take a self-consistent equation of state p ~ {'J — l)e 
with approximate adiabatic index 



7 = 



12p + 5pc^ 
9p + 3pc2 



(5) 



^ After the 2005 outburst this knot is no longer stationary and 
propagates with an apparent velocity of ~ 2.7c (Giovannini et al. 
2010). 



^ For a comprehensive description of the shocked jet zone see 
Kohler et al. (2011). 



© 0000 RAS, MNRAS 000, 000-000 



The Efficiency of Reconfinement Shocks 3 



based on Ryu et al. (2006). 

The local dissipation efficiency is defined as 

/diss /kinj /kiiijS 

Ediss = = 7 , (O) 

/kin,j ,/kin,j 

where /kin ~ (F— l)r/?xpc'^ is the kinetic energy flux density, 
/diss is the dissipated energy flux density and I3± is the di- 
mensionless velocity component perpendicular to the shock 
front. Under the assumption of the cold upstream plasma, 
it can be simplifled to 

. - ^'-^^ (7) 

Ediss - — — ■ [I) 



100% 




3 RESULTS 

In Paper I, we studied the dependence of the average dis- 
sipation efficiency on Fj and Oj for the case of 77 = 0, i.e. 
uniform external pressure. We found that the efficiency de- 
pends sensitively on the product Tj&j. For FjSj < 1, an 
approximate scaling law ediss ~ 6%(rj0j)^ can be used. For 
Fj6j > 1, very high values can be achieved. 

Here, we have additionally calculated the average dis- 
sipation efficiency for different values of the pressure index 
rj. The results are shown in Figure 2. Instead of the opening 
angle 0j, the reconflnement angle Or is used, multiplied by 
Fj, on the horizontal axis. As has been shown by Komis- 
sarov & Falle (1997), the relation between these two angles 
is ©r ~ SOj, where S — 1 — r]/2. For 77 = we have 0r ~ 0j, 
so these results are consistent with the flndings of Paper I. 
We thus generalise the previous result and show that the 
dissipation efficiency is determined by a single parameter in 
the three-dimensional space (rj,Oj,r;). In Figure 2, we also 
plot a slightly different scaling law'', ediss ~ 8%(rj0r)^. A 
noticeable discrepancy for FjOr < 1 can only be seen for the 
case of 77 = 1.5. 

In addition to the average dissipation efficiency, we are 
interested in the longitudinal distribution of the dissipated 
energy. Figure 3 shows the energy flux dissipated per unit 
jet length dz for different external pressure indices 77. The 
proflles of dissipated energy depend very strongly on 77. For 
the flat external pressure distribution (77 = 0), most of 
the energy is dissipated beyond the half of the reconflne- 
ment length. The amount of dissipated energy tends to as 
z — >■ Zr, even though the dissipation efficiency increases with 
z, since the shock becomes less and less oblique and conse- 
quently Fs decreases. But this increase in efficiency is much 
slower than a decrease in the jet cross-section and hence 
a decrease in the jet kinetic energy flux per unit dz. The 
main dissipation region shifts closer to the jet origin with 
increasing 77. The peak of energy dissipation rate Zdiss.max is 
located at ~ O.GS^r for 77 = 0, at ~ O.S^r for 77 = 0.5 and 
at ~ 0.232:r for 77 = 1. For 77 = 1.5, the dissipation profile 
changes to monotonically decreasing with z and the bulk of 
the dissipation takes place very close to the jet origin. 

Figure 4 shows the dependence of Zdiss.max on the jet 

^ The reason for the change of the normalising factor in the 
power-law scaling is that in Paper I the value for rj0j = 1 
has been used, while here we require a good overall match for 
0.1 < TjQr < 1. 



Figure 2. Average dissipation efficiency as a function of the prod- 
uct of the jet Lorentz factor Fj and the reconfinement angle 0r. 
There are 4 families of models for the case of ffat external pressure 
(77 = 0), all plotted with red solid lines, with varying ©j and Fj 
fixed at values 5, 10, 20 and 40. Models for 77 > are calculated 
using Fj = 10. A power-law scaling valid for FjSi- < 1 is shown 
with the dashed black line. 




Figure 3. Energy dissipation rate per unit of jet length dz 
for different external pressure indices rj. The length scale 2 is 
normalised to the reconfinement length z^. The dissipation rate 
profiles are normalised to unity, with -Fdiss = / /dissd^, where 
dS = 27rrad2/ cos Os is the shock front surface area. 



opening angle and the external pressure index. For a given 
77 < 1, the position of the dissipation peak with respect to 
the reconflnement length Zr is relatively stable for FjSj < 
1 and is systematically shifted towards the jet origin for 
FjOj > 1. For 77 = 1.5, the dissipation peak is always lo- 
cated very close to z = 0. The black dashed line shows the 
dependence of the location of the dissipation peak on 77 for 
fixed Fj and 0j. We find that the Zdiss,max/zr ratio decreases 
with 77, falling to the vicinity of for 77 > 1.3. It is remark- 
able that, despite the fact that for 77 > 1 the dissipation is 
concentrated close to the jet origin, the average dissipation 
efficiency scales with the reconfinement angle and not the 
opening angle. And although there are large differences in 
the spatial distributions of the energy dissipation rate, ediss 
scales in the same way for all values of 77. 
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Figure 4. Position z^iss i^^x of the peak of the energy dissipation 
rate along the jet axis, relative to the reconfinement length Zr, as 
a function of the external pressure index 77 and the jet opening 
angle ©j. Solid colour lines show the dependence of 2diss,max/-2r 
on the product of ©j and the jet Lorentz factor Fj (plotted against 
the lower a;-axis) for several values of rj. Dashed black line shows 
the dependence of Zdiss max/^r on jj (plotted against the upper 



X-axis) for Fj = 10 and ©j 



4 DISCUSSION 

The results of this work and Paper I show that relativistic re- 
confinement shocks can be very efficient means of energy dis- 
sipation. Their efficiency depends on a simple combination 
of fundamental parameters of the jet and its environment. 
For many astrophysical jets, their Lorentz factors, opening 
angles and total powers can be measured or significantly con- 
strained. In such cases it is possible to test the hypothesis 
that energy dissipation is dominated by the reconfinement 
shock. 



4.1 GRBs vs AGNs 

Achromatic breaks detected by Swift in some afterglow light 
curves allow one to constrain GRB jet opening angles. In 
several cases it has been found that TjOj ^ 1. It became 
a challenge for numerists studying the initial acceleration 
and colhmation of relativistic jets to reproduce such wide 
jets (Komissarov et al. 2009). The solution was to interrupt 
the coUimation at some point, as would be expected for a 
jet breaking out of its host star (Tchekhovskoy et al. 2010; 
Komissarov et al. 2010). Such a situation is unique for GRB 
jets and allows reconfinement shocks forming at larger dis- 
tances to be very efficient dissipators. 

A scenario for a long GRB involving a very efficient re- 
confinement shock has been investigated numerically by Laz- 
zati et al. (2009). The initial parameters adopted by them 
are Tj = 400 and Gj = 10°, which translates to FjQj = 70. 
Our model predicts in such a case an efficiency of ~ 90% for 
77 = and ~ 82% for 77 = 1. Their numerical result is thus 
consistent with our scaling law. 

Observations of AGN jets indicate that they satisfy the 
relation FjOj < 1 {e.g., Pushkarev et al. 2009), for which 
we predict at most a moderate dissipation efficiency. This 
is consistent with the initial jet coUimation not being in- 
terrupted due to a change in the environment. The value 
of Ediss ~ 8%, corresponding to the case of FjOr ~ 1, is 
in line with estimated radiative efficiencies of the brightest 



blazars. As we show below, the dissipation efficiency of jets 
in low-luminosity AGNs, such as M87, can be much lower. 

4.2 The jet of M87 

Stawarz et al. (2006) provided a thorough review of the 
properties of the M87 jet and its environment. The Lorentz 
factor is estimated at Fj ~ 6 and the viewing angle at 
^obs 20°. The HST-l knot is located at deprojected dis- 
tance of 2r ~ 180 pc. The pressure distribution of the 
host galaxy within zb ~ 230 pc has been estimated as 
Pext(2) = Pb(z/2:b)"'', where pb = 1.5 x 10~® dyncm"^ 
and 7] = 1.2 (hence S = 0.4). The position of the knot and 
the distribution of external pressure can be used to calcu- 
late the jet power, employing a formula derived from the 
analytical model of Komissarov & Falle (1997): 



44 -1 
ergs 



(8) 



where /i = 17/24. This value is a bit higher than the estimate 
10*** ergs^^ obtained from the energetics of the radio lobes 
by BickneU & Begelman (1996). 

This study can be complemented by an estimate of the 
dissipation efficiency. We already know the jet Lorentz fac- 
tor and we need to calculate the reconfinement angle. This is 
complicated by the fact that the jet region immediately up- 
stream the HST-l knot is not visible on any high-resolution 
radio maps. Using the 20 GHz VLBA map from Cheung 
et al. (2007), we measure the projected aspect ratio of the jet 
section up to HST-l: (2rm/zr)proj ~ 0.026. This corresponds 
to the actual aspect ratio of (rm/zr) ~ ('■m/zr)proj sin Sobs ~ 
0.0045. Using the analytical model of Komissarov & Falle 
(1997), it is possible to express the reconfinement angle as 



(9) 



0.83°. Since 



For rj — 1.2, we obtain the value Oi- 
FjOr ~ 0.087 <S 1, the dissipation efficiency can be 
calculated from the approximate scaling relation: ediss ~ 
6 X 10" 



Multiplying it by the jet power estimated in 
Equation 8, we obtain the energy dissipation rate Ldiss = 
Ediss^j ~ 3.3 X lO*'^ ergs"'^. If all of the dissipated en- 
ergy were radiated away, the observed luminosity would 
be Lobs ~ (©?/Fj)Ldiss - 6 X 10*' ergs'^ Here, Vj = 
[Fj(l — /3j cos Sobs)] ~^ ~ 2.3 is the jet Doppler factor and 
the relativistic boost factor (Vf /T^) is used in the form 
that is valid for a stationary emitting region, rather than 
a co-moving one (see Sikora et al. 1997). At the distance 
of dL ~ 16 Mpc, the observed bolometric flux would be 
/obs = Lobs/(47rd^) ~ 2.1 X 10"" ergs-'cm'^ This value 
is a bit higher than the actually observed broad-band flux 
of M87 (see Figure 4 in Abdo et al. 2009), which is domi- 
nated by the non-thermal emission from the inner jet. In- 
terestingly, the observed flux could be matched if the lower 
estimate for the jet power by BickneU & Begelman (1996) 
is used instead of the result of Equation 8. This indicates 
that the simple analytic model of relativistic reconfinement 
shocks overestimates the jet power, but predicts a correct 
dissipation efficiency. 

If the radio emission from the inner jet of M87 is pro- 
duced at the reconfinement shock, its spatial distribution 
should be related to the distribution of the dissipation rate. 
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The radio map from Cheung et ah (2007) shows that emis- 
sion peaks close to the galactic nucleus (see also Hada et al. 
2011) and decreases monotonically with the distance. As 
shown in Figures 3 and 4, the dissipation rate behaves in a 
similar manner for the external pressure index r/ > 1.3. The 
actual index inferred for M87, rj = 1.2, is close to this range. 
Also, the radio map shows an edge-brightened jet, which is a 
natural consequence of dissipation at reconfinement shocks 
(Nalewajko 2009) 



5 CONCLUSIONS 

This work generalises the findings of Nalewajko & Sikora 
(2009, Paper I) on the dissipation eSiciency of relativistic 
reconfinement shocks and sets these studies in a broader 
astrophysical context. 

We find that the average dissipation efficiency depends 
on the product of the jet Lorentz factor Fj and the reconfine- 
ment angle Qr, which is equal to the opening angle Oj for a 
flat distribution of external pressure (77 = 0). For FjSr < 1, 
an approximate scaling law eaiss ~ 8%(rj©r)^ can be used. 
This moderate-efficiency regime can be applied to the jets of 
AGNs, while the high-efficiency regime (TjOr ^ 1) is char- 
acteristic for GRBs. The differences in radiative efficiency 
between these sources may be related to the different cir- 
cumstances of the initial jet collimation process. A similar 
idea for the unification of relativistic jets between GRBs 
and AGNs has been recently formulated by Nemmen et al. 
(2011). 

Our results have been applied to the jet of radio galaxy 
M87, hosting a peculiar knot HST-1. Emission from the in- 
ner jet of M87 is consistent with dissipation at a reconfine- 
ment shock extending upstream from HST-1 in two aspects: 

- the broad-band luminosity of M87 is consistent with 
the product of the dissipation efficiency ediss predicted by 
our model and the independently estimated jet power; 

- radio emission peaking close to the nucleus resembles 
the dissipation profile of reconfinement shocks for external 
pressure index 77 = 1.2. 

Our results suggest that reconfinement shocks may be a 
dominant dissipation mechanism in astrophysical relativistic 
jets. As such, they deserve more attention and more detailed 
investigations. 
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